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Abstract 

Time-dependent thermal X-ray spectra are calculated from physically plausible con¬ 
ditions around GRB. It is shown that account for time-dependent ionization processes 
strongly affects the observed spectra of hot rarefied plasma. These calculations may pro¬ 
vide an alternative explanation to the observed X-ray lines of early GRBs afterglows (such 
as GRB 011211). Our technique will allow one to obtain independent constraints on the 
GRB collimation angle and on the dumpiness of circumstellar matter. 


1 Introduction 

X-ray afterglows of GRBs were first discovered by Beppo-Sax satellite (Costa et al. 1997) 
and their spectral observations are now used for plasma diagnostic of the plasma surrounding 
GRB cites (Lazzati et al. 1999; Piro et al. 2000). The most exciting observations are now 
available from Chandra and XMM-Newton X-ray satellites. Recently, Reeves et al. (2002a,b) 
reported on XMM-Newton observations of the early afterglow spectrum of GRB 011211. These 
observations, taken ~ 11 hours after the GRB, revealed the presence of emission lines which 
were fading over a time scale shorter than 10^ s. These lines were interpreted as blue-shifted 
(with a velocity of ~ 0.1c in the GRB comoving frame) emission from ions of a-elements 
(Mg XI, Si XIV, S XVI, Ar XVIII, Ca XX) arising in a dense (ue ~ lO^^cm”^), optically-thin 
thermal plasma shell heated by the GRB. The model requires a factor of 10 overabundance of 
a-metals in the plasma. No iron lines were seen, however. 

If this interpretation is correct (see the criticism in Borozdin and Trudolyubov 2002; Rut¬ 
ledge and Sako 2002), this could strongly point to a supernova (SN) explosion preceding the 
GRB (the supranova model; Vietri and Stella 1998). The apparent absence of iron lines is ex¬ 
plained by Reeves et al. by possible peculiarities of the SN ejecta. It is worth noting the claims 
(Watson et al. 2002) that three of four GRB X-ray afterglows observed by XMM-Newton are 
best-fitted by thermal plasma model. 
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The emission line origin in GRB 011211 was critically considered in several papers (Kumar 
and Narayan 2002, Lazzati 2002). Our study (which we have reported at this meeting before 
the latter paper was posted in astro-ph e-print archive) re-examines the X-ray emission line 
origin in thermal plasma model. In contrast to other studies, we investigate time-dependent 
ionization effects in thermal plasma. We show inconsistency of the model proposed by Reeves 
et ah (2002a) for GRB 011211 and propose another interpretation of what was observed. 
Specihcally, in our model the observed non-stationarity of emission lines is directly related to 
time-dependent character of collisional ionization in thermal optically thin plasma with clumps 
of lower density (ue ~ 10^^ — 10^^ cm^). Such a plasma can be found around the GRB or just 
spuriously occur near our line of sight. 

However, in our model it is difficult to explain the observed blueshift of lines by a supernova 
preceding the GRB for only a few days. In the thermal model, no individual lines can be 
reproduced from fresh hydrogen-helium SN ejecta even enriched in metals (there is no space 
near the GRB source to get the needed emission measure). Purely metal ejecta are much more 
efficient X-ray emitters, they can be placed closer to the GRB, but then, as our computations 
show, one would rather observe the metal photorecombination edges, not lines, from a recent 
(a few days old) SN ejecta. The ten-times metal overabundance is also unnatural for fresh 
ejecta, because modern paradigm of massive star SN explosions assumes purely metal ejecta 
(for SNe Ic). If real, the blueshifted lines should be emitted at much larger distance from 
the GRB engine, so they could originate in the material expelled by a supernova much earlier, 
maybe some years before the GRB, and this material can be mixed with the H-rich circumstellar 
matter. 


2 Time-dependent emission lines from thermal plasma 


The thermal emission model assumes collisional ionization of atoms and is applied to rarehed 
plasmas heated by an external source (the solar corona provides an example). Photons orig¬ 
inating in atomic transitions, in radiative recombination, etc. leave the plasma freely. The 
plasma cooling time depends on density and chemical composition; for example, for a fully 
ionized plasma composed of ions with the charge Z, 
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The numerical factor in this formula is valid with the 20%-accuracy for a purely hydrogen 
plasma. Glearly, this formula gives an upper limit for the cooling time in real plasmas. 

For densities Ug ~ lO^^cm”^, advocated by Reeves et ah, the plasma cools faster than 
in 1 second and a stationary thermal model is appropriate. At lower densities, for example. 
He ~ 10^^ — 10^^ cm“^, the cooling time becomes much longer than the heating time th ~ 90 
s (the GRB 011211 duration in its rest frame) and time-dependent effects must be taken into 
account. The relevant parameter here is nt. For large nt, plasma relaxes to a stationary state. 

In the case of GRB X-ray afterglows, the non-stationarity has two faces. On one hand, the 
real X-ray spectrum could change over the integration time in the detector (~ 5000 s for GRB 
012111) due to temperature decrease even if the plasma density is high and the parameter nt 
formally allows for stationary ionization. On the other hand, when nt is small, relaxation effects 
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could cause ions to appear which were absent in the stationary case. We stress here that there 
is strong dependence of plasma cooling time on the metallicity. The real cooling time depends 
on ionization states and abundances of different ions so time-dependent effects play signihcant 
role in low-density plasma. 
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Figure 1: Relative number densities of ions as a function of temperature for stationary case 
(left panel) and as a function of time for time-dependent model (right panel). 

This latter face of non-stationarity is illustrated in FigQ which shows relative number den¬ 
sities of some ions as a function of temperature (for stationary case) and as a function of the 
parameter nt for time-dependent ionization. In the latter case, an initially cold plasma (10® 
K) with density n = 10^^ cm“® is instantly heated up to a temperature of 10® K. This hgure 
indicates that during thermal relaxation of plasma, ionization state at a given time can be close 
to some stationary ionization state at different hxed temperature. 

3 X-ray afterglow of GRB 011211 

X-ray afterglow of GRB 011211 was observed by XMM-Newton ~ 11 h after the GRB (Reeves 
et ah 2002a,b). The emission lines faded away over a time interval of 10^ s. The volume 
emission measure as inferred by Reeves et ah from the average X-ray luminosity 7 x 10^® erg/s 
for plasma with normal cosmic abundance is EM = rigR = 10®®cm“^. As we discussed above, 
such an estimate is only the upper limit, since real plasma can have higher volume emissivity 
and smaller cooling time. 

3.1 Thermal model by Reeves et al. 

According to the model by Reeves et al. (2002a), the observed X-ray spectrum originates in a 
dense Ug = 10^® cm“® shell of radius R = 10^® cm and thickness h = 10® cm expanding with a 
velocity of 0.1 c from GRB, as the mean redshift of the identihed lines {z = 1.88) implies (the 
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GRB host galaxy redshift is 2 : = 2.14). As such a dense plasma cools down in 1 s, the observed 
fading time of X-ray lines (tj = 10^ s) is explained by the light retardation from different parts 
of instantly illuminated shell (this geometrical effect was invoked by Lazzati et al. 1999 to 
explain the X-ray afterglow evolution of other GRBs). 

This model has obvious caveats which make its applying to GRB 011211 questionable. 

• Short plasma cooling time: 

If cooling time is as short as 1 s, an overwhelming part of the observed X-ray emission 
must be generated during illumination of the plasma by gamma-ray radiation from the 
GRB {th = 90 s in the GRB rest frame). The time evolution is then entirely due to 
the retardation, i.e. to the geometrical effect. This implies that the observed X-ray 
photons were actually generated simultaneously with the powerful illumination by hard 
photons fr-om the GRB. Then, in order to explain the corresponding X-ray spectrum 
the photoionization (nebular) approximation but not the thermal (coronal) one should 
be used. However, the photoionization model is claimed to be inadequate (Reeves et al. 
2002a). In order to try to explain these observations within the frame of a thermal plasma 
model, we must assume much longer cooling time tc> thi hence lower plasma densities. 

• Large emission measure: 

The EM in this model is determined from the mean X-ray continuum luminosity and with 
the assumed density Ue = 10^® cm“^ relates to the layer volume V = R^h = (10^^)^ x 10® = 
10^® cm^ which we observe in 1 s. As the X-ray luminosity is registered over 10^ s, there 
must be 10^ such volumes! This contradicts the assumption that one and the same shell 
of radius R = 10^^ cm and thickness h = 10® cm is emmitting. So in order to save the 
model, on has to assume the existence of clumps with size I = 10^^ cm 3> h, which are 
optically thick at such densities. 

• Apparent absence of iron lines: 

The explanation was that iron ^®Fe does not have enough time to be produced from the 
radioactive ^®Ni expelled by SN. Under special assumptions about SN II mechanism one 
can obtain the dehcit of iron isotopes (Rausher et al. 2001). However, if the plasma 
shell in this model is indeed expelled by SN, such an SN is hardly to be typical - in 
particular, it can be asymmetrical. So the initial amounts of iron, nickel, and cobalt 
can be comparable. Asymmetric explosions of SN that may be connected to GRBs give 
different results (Pruet et al. 2002). Some model calculations of hypernovae by Umeda 
et al. (2002) yield up to 40% of prompt iron isotopes. So the iron yields in SN ejecta are 
strongly model-dependent, and the apparent absence of iron in the X-ray spectrum can 
not provide a reliable argument in favor of a supernova associated with the GRB. 

• High density and high velocity of the shell: 

Physical characteristics of the SN ejecta as a dense rapidly expanding shell also meet 
difficulties. Similar dense shells are produced in radiation hydrodynamics calculations of 
SN (Falk and Arnett 1977; Ghevalier and Klein 1979; Blinnikov and Bartunov 1993), but 
their properties are quite different. The total energy of the shell in the model by Reeves 
et al. approaches the typical SN explosion energies, while the dense shells arising in hydro 
calculations either have smaller masses or lower velocities, of order of several 1000 km/s 
(Ghugai et al. 2002) and carry only a small fraction of the SN explosion energy. 
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3.2 Thermal model with time-dependent ionization 

In our model, the observed X-ray fading time is identified with thermal relaxation time of a 
plasma heated up by GRB. 

3.2.1 Emitting region geometry 

Let the GRB occur inside (or nearby) a star forming region of complex clumpy structure. This 
assumption is quite natural and has solid observational support (van Paradijs et ah 2000). In 
order for plasma cooling time (HD to be 10^ s, the gas density must be Ue ~ 10^^ — 10^^ cm 
To obtain the observed EM (n^ V < 10®® cm“®, the volume of the emitting region must be 
V < 10*^^ — 10^® cm^. At the same time, it must be optically thin for X- and gamma-ray 
photons, nlaT ~ 1, where ax ~ 10“^'^ cm^ is Thompson electron cross-section. It is hard to 
meet both conditions in a homogeneous medium, but a more realistic rarehed plasma with 
clumpy structures can do the job. 

Indeed, let us consider ~ 10® gas clouds with density ~ 10^^ cm“® each and a size ~ 
10^® cm. The total volume of the dense component in such clouds amounts to the required 
value 10^® cm® and gives the necessary EM = 10®® cm“® (we recall that these are upper limits; 
for example, using realistic plasma cooling function with normal cosmic abundances yields 
EM ~ 4 X 10®® cm“® and solar plasma without H and He yields EM ~ 10®®cm“®, see also 
FigEJ. The total mass of the emitting clumps is less than ~ O.SMq. Glumps of such densities 
and size are indeed observed as maser condensations in stellar winds from young stars and late- 
type stars (e.g. Bains et ah 2002, Richards et ah 1999, Lekht et ah 2001). The total volume 
occupied by this inhomogeneous medium is determined by energy balance and the width of a 
gamma-ray beam. 




Figure 2: Plasma cooling function (erg cm® s“®). Left panel; solar composition, then the 
emission measure for T = 10® K is EM = ri^V = Lxl£~ — 3.7 x 10®® cm®; Right panel: 
metal-rich medium (without H and He), EM = ri^V = Lx/€~ ~ 1.4 x 10®® cm® 

Heating matter up to 10® K would require energy deposit cq = lOkeV per nucleon, i.e. 
the energy Ep = 10^^cm“® x 10^® cm® x 10 keV = 10^® erg is necessary. The total isotropic 
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energy of GRB 011211 is = 5 x 10^^ ergs. The effective cross-section of gamma-ray energy 
deposition per nucleon with good accuracy is ~ O.lcrr (Ambwani and Sutherland 1988). 
Considering the fluence from the source at a distance d to be F = the heating energy 

per nucleon is Cq = so 

d = ~ 2 X 10'" cm , 

V dTreo 

i.e. the GRB energy is sufficient to heat up the surrounding matter up to distances d ~ 0.1 
pc. Note that F plays more important role than E^] it is measured by our detectors (at larger 
distances of course), and no assumptions on gamma-ray emission beaming should be made to 
estimate d. 

The gamma-ray beam opening angle 6 determines the total volume of the illuminated mat¬ 
ter. The simple geometry leads to the estimate 

d(l — cosd) = ctj/(l -I- x) ~ lO'^ cm 

Assuming not very large angle (which appears natural in view of the large value of isotropic 
E-y) we have d(l — cosd) ~ d6‘^/2, i.e. the illuminating gamma-ray beam size is R cx 6d = 
(2dctj/(l -1- x))'A = (2 X 10'^ cm x d)'A. Por example, R ~ 10'® cm for d = 2 x 10'" cm and 
6 cx Rld = 0.05, i.e. 3° (the total beam angle is twice as large, 29 = 6° ). 

The total volume of the medium illuminated by GRB is quite large, V R^d = 3x 10^9 
cm“^. We stress here that though gamma-ray photons from GRB fly this distance d in several 
months, they reach the detector simultaneously with X-ray photons they created passing the 
clouds, so X-ray emission comes not from a dense shell 10'"' cm in thickness but from an 
extended clumpy region with the size of order of d. The hlling factor of this medium can be 
10“^ or larger if the condensations are located in a relatively thin layer and not across the entire 
distance d. This simple estimate shows that the physically acceptable structure of a clumpy 
ISM can actually be the source of the observed X-ray emission. 

The necessity of clumpy structure of emitting plasma was also obtained by Lazzati (2002) 
from independent analysis of stationary thermal X-ray afterglows from GRBs. We stress again 
here that in contrast to Lazzati, we do not require the shell around the GRB, for us it suffices 
that the medium with the required physical parameters lies in the cone of the gamma-ray beam. 
This plasma can be not genetically connected with GRB. In our model, orphan X-ray afterglows 
can well show transient emission features like those observed in the case of GRB 011211. 

3.2.2 Spectral modelling 

To model X-ray spectra from thermal plasma, including the case of GRB 011211, we employ a 
numerical code originally worked out by P. Lundqvist for calculations of collisional ionization 
of a stationary plasma, which takes into account the following basic processes: ionization by 
electron impact, autoionization, photorecombination, dielectronic recombination, and charge 
transfer. The algorithm outputs bremsstrahlung (free-free), free-bound, bound-bound, and 
two-photon emission. The code has been upgraded by S.I. Blinnikov and E.I. Sorokina to a 
time-dependent version {hi = f{T{t),n{t),Z,...) for all ionic species). Input parameters are 
initial temperature, density, and abundances of elements. 

To calculate time-dependent spectra, the temperature dependence on time T{t) is needed. 
The plasma cooling function (erg cm^ s“') is almost independent of temperature for 2 x 10®K 
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- 10®K (see Figl2I). Then the temperature change with time can be found from the differential 
equation 

dT _ e~{To)nl 
dt Cvp 

where Cy = A+X^.) is heat capacity per unit mass at hxed volume, TZ is the universal 

gas constant, A is the mean mass of ions in atomic units, Xf, is the number of free electrons 
per nucleon, Tq = 10®K is the initial temperature. Then in the above temperature range the 
temperature decreases linearly with time 

T{t) = To-e-{To)nlt/{Cvp) (2) 
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Figure 3: Relative number densities of ions as a function of the parameter nt for linear plasma 
cooling law Q. Left panel: solar composition; right panel: metal-rich plasma (without H 
and He) 

The time evolution of ion number densities (FigO)) is calculated numerically and determines 
the appearance of the corresponding emission X-ray lines (FigE)). These hgures show strong 
time evolution of thermal spectra from optically thin cooling plasma initially heated by GRB 
up to 10-keV temperatures. For normal cosmic abundances, emission lines from moderately 
ionized iron atoms immediately appear. As plasma cools down, lines of H-like and He-like iron 
emerge (Figl^. If metal enriched (and, correspondingly, with deficit of H and He) SN ejecta is 
observed, emission lines from metals become more prominent against the continuum. At late 
cooling stages, photorecombination edges of metals prevail (FigH)). 


4 Discussion 


Analysis of our calculations of emission spectra from thermal optically thin plasma suggests 
that interpretation of early X-ray GRB afterglows in term of thermal model requires (1) a 
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Figure 4: Time evolution of X-ray spectrum (in the observer frame at the GRB redshift z = 
2.14) for different time moments. Left panel: solar composition; right panel: metal-rich 
plasma (without H and He) 

moderate density of the emitting plasma {ue ~ 10^^ —10^^ cm“^), and (2) a clumpy structure of 
the circumstellar (interstellar) medium. In our model, X-ray line emission forms not necessarily 
in matter genetically related to gamma-ray burst source, it is sufficient for it to he within the 
gamma-ray illumination beam at distances < 0.1 pc from the GRB site. Specihcally, for the 
observed parameters of GRB 011211, Lx(0.6 — 30) keV 7 x lO'^^ erg/s (in the source rest- 
frame), a total of ~ O. 3 M 0 of matter comprised in clouds ~ 1 AU across inside a volume of 
~ 10^® cm^ illuminated by gamma-ray beam with the full width 29 ~ 6° would suffice. Note 
that such narrow angles are generally obtained in the standard model interpretation of time 
decay of GRB afterglows (Frail et ah 2001). The parameters of individual clouds are also 
astrophysically feasible and correspond to dense maser condensations in stellar winds from 
AGB stars (Bains et ah 2002, Richards et ah 1999, Lekht et ah 2001). 

The proposed model interprets the observed temporal evolution of X-ray spectrum in terms 
of thermal relaxation time of ions and plasma cooling, and does not appeal to purely geometrical 
arguments. However, the comparison of the calculated spectra with those observed by Reeves 
et ah does require enhanced a-metal abundances. The apparent absence of iron lines either 
requires its underabundance, which is less probable according to current SN nucleosynthesis, 
or can be explained by their rapid variability or washing out against a high level of nonthermal 
power-law continuum (which is not excluded by spectral htting). 

Thermal X-ray spectra of a metal-rich shell ejected by an SN preceding GRB for a few 
days would rather show photo-recombination edges and not stationary emission lines with time 
evolution due to geometrical factor. 

In our model, it appears more natural to consider that no blue-shifted X-ray lines were 
observed, if the clumpy medium is a circumstellar envelope produced by slow superwind from 
the GRB progenitor. The case is even stronger if what we observe is a nearby clumpy cloud in 
a star-forming region. If the signihcant blue-shift of the X-ray lines will be conhrmed, it will 
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provide a strong evidence that the GRB is preceded by a supernova. However, if these rapidly 
moving clouds are hydrogen rich, they can hardly be located at such small distances as inferred 
by Reeves et al. (2002a) and the preceding supernova must occur several months or even years 
before the GRB event (the supranova scenario, Vietri and Stella 1998). 

Apparently, the adequate explanation of the observed spectrum of X-ray afterglow of GRB 
011211 should be sought for within the frames of a combined model of photoionization heating 
of an inhomogeneous medium around the GRB progenitor and time-dependent collisional ion¬ 
ization in more dense condensations with an optical thickness for scattering about 1. Glearly, 
the calculation of realistic X-ray spectra from such media requires appropriate account for 
radiative transfer effects, which are beyond the scope of the present work. 

5 Conclusion 

We have calculated thermal X-ray emission spectra of optically thin thermal plasma using a 
time-dependent treatment of collisional ionization. We show that an instant heating of this 
plasma up to the 10-keV temperatures produces time-dependent X-ray emission spectra in 
which lines of a-element metal ions and iron at different ionization states are apparent. We 
propose a model of thermal X-ray emission for the observed early X-ray afterglow of GRB 
011211 , in which a narrow gamma-ray beam passes through clumpy interstellar medium, not 
obligatorily genetically related to GRB. The non-stationarity of the X-ray spectrum is related 
not to a geometrical factor (time retardation from different parts of the shell), but to a physical 
non-stationarity of the collisional ionization and thermal relaxation of plasma with density of 
order 10^^ — 10^^ cm“^. Parameters of such a medium resemble maser condensations in active 
star-forming regions. We conclude that the account for time-dependent effects in plasma is 
unavoidable in spectral X-ray diagnostics of the medium around GRB sites within the frame 
of thermal plasma models. 
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through RFBR grants 02-02-27445z and 01-15-99310m. S.I. Blinnikov and E.I. Sorokina are 
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